

if("`c(os)'"=="Windows"){
	local tableSaveDir = "Tables\"
	local plotSaveDir = "Plots\"
	local pythonSaveDir = "PythonScripts\"
}
else{
	local tableSaveDir = "Tables/"
	local plotSaveDir = "Plots/"
	local pythonSaveDir = "PythonScripts/"
}




local numIter = 100

matrix betas = [0]

forvalues iter = 1(1)`numIter'{
	clear
	import excel "Data/SimData.xlsx", sheet("Travel-`iter'") first
	
	reghdfe Fleet TreatInd, absorb(State Time)
	local pi = e(b)[1,1]
	reghdfe Use TreatInd, absorb(State Time)
	local delta = e(b)[1,1]
	local beta = `delta'/`pi'
	
	di "`delta'"
	di "`pi'"
	di "`beta'"
	
	matrix betas = [betas \ `beta']
	
}


matrix betas = betas[2..., 1]

matlist betas


matrix gammas = [0]

forvalues iter = 1(1)`numIter'{
	clear
	import excel "Data/SimData.xlsx", sheet("Gas-`iter'") first
	
	reghdfe Fleet TreatInd, absorb(State Time)
	local pi = e(b)[1,1]
	reghdfe Use TreatInd, absorb(State Time)
	local delta = e(b)[1,1]
	local beta = `delta'/`pi'
	
	di "`delta'"
	di "`pi'"
	di "`beta'"
	
	matrix gammas = [gammas \ `beta']
	
}


matrix gammas = gammas[2..., 1]

matlist gammas


svmat gammas
svmat betas

histogram gammas1, title("Simulated Gammas") xtitle("")
graph export "Plots/AdditionalSimulations/Gammas.png", replace
histogram betas1, title("Simulated Betas") xtitle("")
graph export "Plots/AdditionalSimulations/Betas.png", replace
